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Non-perturbative approaches to quantum gravity call for a deep understanding of the emergence 
of geometry and locality from the quantum state of the gravitational field. Without background 
geometry, the notion of distance should entirely emerge from the correlations between the gravity 
fluctuations. In the context of loop quantum gravity, quantum states of geometry are defined as 
spin networks. These are graphs decorated with spin and intertwiners, which represent quantized 
excitations of areas and volumes of the space geometry. Here, we develop the condensed matter 
point of view on extracting the physical and geometrical information out of spin network states: we 
introduce new Ising spin network states, both in 2d on a square lattice and in 3d on a hexagonal 
lattice, whose correlations map onto the usual Ising model in statistical physics. We construct these 
states from the basic holonomy operators of loop gravity and derive a set of local Hamiltonian 
constraints which entirely characterize our states. We discuss their phase diagram and show how 
the distance can be reconstructed from the correlations in the various phases. Finally, we propose 
generalizations of these Ising states, which open the perspective to study the coarse graining and 
dynamics of spin network states using well-known condensed matter techniques and results. 


I. INTRODUCTION 

The search for a quantum theory of gravity remains 
one of the major challenge of theoretical physics. Even if 
a complete theory of quantum gravity is still elusive, nu¬ 
merous approaches have shed lights on the problem itself 
and gave insights on the possible structure of the final 
theory. In general relativity, gravity does not evolve in 
a pre-assumed background geometry but moulds the ge¬ 
ometry itself: the gravitational field is encoded directly 
as a curved and dynamical geometry. This background 
independence, conjugated to the major difficulty of con¬ 
structing local diffeomorphism-invariant observables (see 
e.g. [IH3]), leads to the problem of reconstructing the ge¬ 
ometry in quantum gravity from the information stored 
in the quantum state of the gravitational field. 

The line of research we pursue here is at the inter¬ 
play between condensed matter and quantum informa¬ 
tion on one side and quantum gravity on the other: the 
goal is to understand how the distance can be recovered 
from correlation and entanglement between sub-systems 
of the quantum gravity state. The issue of the emer¬ 
gence of distance is of course one aspect of a more gen¬ 
eral program to understand the emergence of geometry 
at a semi-classical level from a purely background inde¬ 
pendent quantum theory of gravity. To reconstruct a 
notion of locality and the geometry, the distance seems 
nevertheless to be the first natural quantity to define. 
Thus focusing on the distance, one could extract it from 
the 2-point correlation function assuming that the inverse 
square law holds. This perspective is quite similar to the 
point of view of spectral geometry and non-commutative 
geometry, in which distances, and more generally the ge- 
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ometry, emerge from the spectral analysis of the 2-points 
function, given respectively by the Laplacian and Dirac 
operators mg. Like in condensed matter physics, we 
expect the behavior of the correlations to depend on the 
considered phase. While the typical behavior of correla¬ 
tions is an exponential decay, the interesting regime for 
gravity, namely an algebraic decay, would be attained in 
critical regimes or algebraic phases of some specific model 
like the XY model. In this paper, we would like to pro¬ 
pose a simple framework to test these ideas in the con¬ 
text of loop quantum gravity and we will introduce toy 
quantum states, which we dub Ising spin network states, 
whose correlations and phase diagram can be mapped 
onto the standard Ising model in statistical physics. This 
will allow us to discuss the relation between correlation 
and geometry in a well-under-control environment. 

Loop quantum gravity is one of the most devel¬ 
oped approach to quantum gravity and proposes a non- 
pertubative canonical quantization of general relativity 
(for textbooks, see [IHH]). The physical picture that 
emerges from this theory describes quantum space-time 
as a discrete geometry where areas and volumes opera¬ 
tors have a discrete spectrum. The formalism is based 
on a 3 -|- 1 formulation of general relativity in terms of 
the Ashtekar-Barbero connection and a densitized triad, 
thus describing gravity as SU{2) gauge field theory, in 
fact a topological constrained field theory. In the quan¬ 
tum theory, spin network states form a basis of the wave 
functions that describes the quantum state of 3d space at 
the kinematical level. Physical spin network states will be 
the solution of the Hamiltonian constraints, which imple¬ 
ment the Einstein equations and generate the dynamics. 

The core of our investigation will be on correlations and 
entanglement entropy on spin network states. In the long 
run, the aim is to analyze the typical -possibly universal- 
structure of spin network states with well-behaved al¬ 
gebraic correlations and satisfying an area law entropy. 
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Correlations and especially entropy are of special impor¬ 
tance for the understanding of black holes dynamics. In¬ 
deed, understanding the microscopic origin of black holes 
entropy is one the major test of any attempt to quantify 
gravity and entanglement between the horizon and its en¬ 
vironment degrees of freedom appears crucial (for recent 
developments in loop quantum gravity, see HUH]). The 
usual spin network basis states are eigenstates of the area 
and volume operators, but they carry trivial correlations 
because of their factorized structure with no entangle¬ 
ment at all between intertwiners degrees of freedom. We 
therefore look for more involved spin network states, cor¬ 
relating the spin and intertwiner degrees of freedom living 
on the whole graph. Having in mind statistical and con¬ 
densed matter physics where studies of correlations are 
fairly advanced both numerically and analytically HU, 
the method we propose is to construct pure spin net¬ 
work states from well-known condensed matter models 
and possibly highlights some relevant lessons for loop 
quantum gravity. A similar endeavor is currently being 
investigated by Bianchi and collaborators, similarly using 
condensed matter methods, who introduce squeezed spin 
network states by a Bogoliubov transformation and focus 
on the evaluation of their entanglement entropy HU. 

Ideally, testing the relation between correlations and 
distance requires physical states solutions of the Hamilto¬ 
nian constraints of loop quantum gravity. However, these 
constraints are rather complicated and, despite several 
proposals, no consensus has been reached on the struc¬ 
ture of physical states solving those constraints at the 
quantum level (see e.g. m for a review of the recent re¬ 
search efforts on the quantum dynamics of loop gravity). 
To circumvent this issue, we propose a different strategy: 
build test spin network states on regular lattices, with 
an a priori notion of distance, and compare that natural 
lattice distance with the reconstructed distance emerg¬ 
ing from the correlation carried by the quantum states. 
Working in such a controlled setting should allow us to 
thoroughly illustrate the link between distance and corre¬ 
lation, test the viability of the proposal of reconstructing 
the distance entirely from correlations and understand 
the structures of well-behaved states supporting this pro¬ 
posal. 

A spin network state is defined on a graph, dressed 
with spins on the edge and intertwiners at the vertices. 
A spin on an edge e is a half-integer je G N/2 giving an 
irreducible representation of SU(2) while an intertwiner 
at a vertex v is an invariant tensor, or singlet state, be¬ 
tween the representations living on the edges attached 
to that vertex. Spins and intertwiners respectively carry 
the basic quanta of area and volume. We build our spin 
network states based on three clear simplifications: 

1. We work on a fixed graph, discarding graph su¬ 
perposition and graph changing dynamics for now, 
and we will focus on working with a fixed regular 
lattice. 

2. We freeze all the spins on all the graph edges. We 


fix them to their smallest possible value, which 
correspond to the most basic excitation of geom¬ 
etry in loop quantum gravity, thus representing a 
quantum geometry directly at the Planck scale. 

3. We restrict ourselves to 4-valent vertices, which 
represent the basic quanta of volume in loop quan¬ 
tum gravity, dual to quantum tetrahedra. 

Since the spins are frozen, the only degrees of freedom 
left are the intertwiners living at the 4-valent vertices. At 
each vertex, these recouple between 4 spins which cor¬ 
responds to a two-dimensional Hilbert space or two-level 
quantum system with pairs of spins recoupling to a spin 
0 or 1. This is completely equivalent to working with a 
qubit at each vertex. Qubits, as the elementary physical 
systems encoding information, are the basic ingredients 
of quantum information and quantum computing and are 
thus the basic building blocks of condensed matter mod¬ 
els. Therefore these simplifications provide us with the 
perfect setting to map spin network states, describing 
the Planck scale quantum geometry, to qubit-based con¬ 
densed matter models. Such models have been exten¬ 
sively studied in statistical physics and much is known 
on their phase diagrams and correlation functions, and 
we hope to be able to import these results to the context 
of loop quantum gravity. One of the most useful model 
is the Ising model whose relevance goes from modeling 
binary mixture to the magnetism of matter. We thus 
naturally propose to construct and investigate Ising spin 
network states. 

In three space dimensions, the natural 4-valent regular 
lattice is the diamond lattice, which is a honeycomb lat¬ 
tice and can be seen as dual to a tetrahedral discretiza¬ 
tion of space. It is different from the usual cubic lat¬ 
tice used in lattice gauge theory, discretized gravity (a la 
Regge) or numerical general relativity, but we believe it is 
more natural and better suited to the loop quantum grav¬ 
ity context. This is nevertheless a brand new proposal to 
construct spin networks on the regular 3d diamond lat¬ 
tice. We will define 3d Ising spin network states on that 
lattice and analyze their properties. For mathematical 
simplicity and representation purposes, we will first focus 
on their two-dimensional counterpart and define 2d Ising 
spin network states on a square lattice. In fact, many 
obtained results are straightforwardly generalized to 3d. 
The 2d Ising model on a regular square lattice without 
magnetic field has been solved exactly m and exhibits 
a phase transition with algebraically decaying correla¬ 
tions at the critical coupling. The phase diagram can be 
studied using quantum field theory methods. Since away 
from the transition the correlations decay exponentially, 
the relevant scenario for gravity is indeed at the critical 
coupling which also has the major advantage to possess 
a continuum limit. 

Thus this state appears to fall into the category of the 
test states we are looking for. Moreover, this method of 
investigation allows us to understand the typical struc¬ 
ture of such states in terms of the holonomy operators 
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of loop quantum gravity: non-trivial correlations seem 
to require holonomies around arbitrarily extended loops. 
We also show that this state is a unique solution of a 
set of local Hamiltonian constraints illustrating how local 
constraints can give rise to extended holonomies and to 
long range correlations, a common feature of condensed 
matter systems. 

The paper is structured as follows. Section |TT] reviews 
the definition of spin network and analyze the structure 
of 4-valent intertwiners between spins ^ leading to the 
effective two-state systems which we use to define the 
Ising spin network states. Different equivalent definitions 
are given in terms of the high and low temperature ex¬ 
pansions of the Ising model. The loop representation of 
the spin network is then obtained and studied as well 
as the associated density which gives information about 


parallel transport in the classical limit. Section III 


troduces a set of local Hamiltonian constraints for which 
the Ising state is a unique solution and elaborates on 
their usefulness for understanding the coarse-graining of 
spin network m and the dynamic of loop quantum grav¬ 
ity. Section |IV| discusses the phase diagram of our Ising 
states and their continuum limit as well as the distance 
from correlation point of view. Section [V| end this paper 
with discussion and perspective. 


each vertex v lives an SU{2) invariant tensor called inter¬ 
twiner iy G Inv 5 [/( 2 ) ((g)e 3 i>V^=), with is the (2je + l)- 
dimensional Hilbert space for the spin jg representation. 
These spin network basis states, as illustrated on figj^ 
define the basic excitations of the quantum geometry and 
they are provided with a natural interpretation in terms 
of discrete geometry with the spins giving the quanta of 
area and the intertwiners giving the quanta of volume. 
This interpretation is most clear in the twisted geometry 
picture m- 



FIG. 1 . Example of a colored graph F with spins ji.2,3 at the 
source node v where the interwiner belongs. The shaded 
region represent an example of the dual geometry encoded in 
the spin network state. 


II. ISING SPIN NETWORK STATE 
A. Definition 

In loop quantum gravity, kinematical states are cylin¬ 
drical functionals of the Ashtekar-Barbero SU(2) con¬ 
nection that depend only on its holonomies along a set 
of paths which define an oriented graph T. We call E 
the number of edges of te graph T and V its number 
of vertices. These wave-functions of the geometry de¬ 
pend on E group elements in SU(2). We further re¬ 
quire those states to be gauge invariant, which means 
an invariance under SU(2) transformations at each ver¬ 
tex. The Hilbert space of wave-functions on the graph 
r is then L‘^{SU{2)^/SU{2)^, d/r) where the scalar prod¬ 
uct is defined by the SU(2) Haar measure dfi. A basis 
of this Hilbert space is given by spin networks states 
using the Peter-Weyl theorem stating that a basis of 
L^(SU(2),d/r) is given by the Wigner matrices of the 
SU(2) group element in all possible irreducible unitary 
representations. On each link e of the graph lives a 
spin jf. associated to the area degrees of freedom. At 


^ In fact, the true states of geometry in loop quantum gravity are 
defined as a sum over all possible graphs F thanks to a pro¬ 
jective limit llfil of graphs Rkin = lim Hp. The measure on 

r—f-t-00 

this space of generalized connections A is given by the Ashtekar- 
Lewandowski measure So iFkin ~ d/r^^). This 

projective limit allows to make sense rigorously of superposition 
of states living on different graphs. 


More explicitly, the spin network state jT, je, iy), based 
on the graph T and dressed with spins je and intertwiners 
iy, defines a wave function '4’{9e) on the space of discrete 
connections of SU{2)^/SU(2)^ 


G) (1) 

e V 


where the trace contracts the intertwiners with the 
Wigner matrix representing the group element 

ge along the links between vertices. The scalar product 
between those states are obtained using the orthogonality 
relation of the Wigner matrices: 




Those spin network states, being tensor product of in¬ 
tertwiners, have trivial correlations. We would like to 
build quantum superposition of those basis spin network 
states that would cary non-trivial correlation. What we 
will develop below is a method to construct such states 
whose correlation structure maps onto the Ising model. 
We will then be able to translate the known quantities 
and results of statistical physics to the context of gravity. 

As outlined in the introduction, we propose to define 
Ising spin network states by freezing the spins to their 
smallest value, je = ^l‘2, that is the smallest quanta of 
area at the Planck scale, and to focus on 4-valent inter¬ 
twiners. These 4-valent vertices will be organized along 
a regular lattice. We will consider the 3d diamond lattice 
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and the 2d square lattice. We will focus on the square 
lattice as first step of analyzing the 3d case. Since most 
of the results obtained in 2d are straightforwardly gen¬ 
eralized to 3d as we will see, we postpone the discussion 
of the 3d Ising spin network to the later section |IVD[ 
Those regular lattices carry a natural a priori notion of 
distance but nothing says it corresponds to a physical 
distance. We will strive to compare this lattice distance 
with the emergent distance reconstructed from the phys¬ 
ical correlation carried by the quantum state. 

1/2 1/2 1/2 


1/2 - 1 )- 1/2 -> 


1/2 1/2 1/2 

FIG. 2. Graphic representation of the 4-valent spin 1/2 and 
one of its possible decomposition into a spin 0 and 1 states. 
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4 4 4 


s channel t channel u channel 

FIG. 3. Graphic representation of the different decomposi¬ 
tions of the 4-valent spin 1/2 intertwiner. A spin 0 or 1 state 
can be associated to each of these decompositions. 


quanta defined by the intertwiners. This will provide the 
geometrical interpretation of our spin network states as 
excitations of volumes located at each lattice node. For 
a 4-valent vertex, this operator is defined as: 


Uv 



{SirGh^f (-i[Ji./ 2 , J 1 .J 3 ]) 


(5) 


In this setting, the space of 4-valent interwiners be¬ 
tween four spins 1/2 is two dimensional: it can be de¬ 
composed into spin 0 and spin 1 states by combining the 
spins by pairs, as on figure Different such decomposi¬ 
tions exist and figure gives a graphical representation 
of them. There are three such decompositions, depend¬ 
ing on which spins are paired together, which we dub s, 
t and u channels as in particle physics and quantum field 
theory. Writing V for the 2-dimensional Hilbert space of 
a spin 1/2, we denote its basis states (|t), |4,)). We call 
I Os) and |ls) the spin 0 and 1 states in the s channel. 
They can be explicitly written in terms of the up and 
down states of the four spins: 

|0s) = I (inn) + imt) - itut) - inn)) 

|i«) = ^ (ittu) + intt) 

_ inn) + inn) + inn) + inn) ^ 


Those two states obviously form a basis of the intertwiner 
space. We give can the transformation matrices between 
this basis and the two other channels: 


|0t) 

| 1 *) 

| 0 „) 

in) 


1/1 V3\ /|0s) 

2 -ij (vin) 

1 / -1 V 3 \ /|0s) 

2 1-73 -ij I in) 


(4) 


The intertwiner basis that actually interests us is not at¬ 
tached to one of these channels but is defined in terms of 
the square volume operator Uy of loop quantum gravity. 
Since the spins, and thus area quanta, are fixed, the only 
freedom left in the spin network states are the volume 


where Ji are the spin operators acting on the i link. For 
clarity, we will for the remaining of this paper use the 
natural quantum gravity units for which G = c = h = 1. 
This operator is Hermitian but not positive; it also reg¬ 
isters the space orientation (or more precisely the orien¬ 
tation of the bivector space over which is itself iso¬ 
morphic to K^). The volume itself can then obtained by 
taking the squareroot of the absolute value of U. Geo¬ 
metrically, 4-valent intertwiners are interpreted as repre¬ 
senting quantum tetrahedron, which become the building 
block of the quantum geometry in loop quantum gravity 
and spinfoam models HHiiis]. Using the change of basis 
relations Uy takes the following form in the s channel 
basis: 




( 6 ) 


The eigenstates Im^,^) can be obtained directly has 
|u^_^) = W (|0g) ipi|l^)) with associated eigenvalues in 

Planck units . This (square) volume is the small¬ 

est non trivial possible value of a chunk of space. The 
only freedom is the orientation of the chunk of volume 
corresponding to the intertwiner state. 

We consider these two oriented volume states \u^^i) as 
the two levels of an effective qubit. Let us point out 
that these two level states are a priori only defined up 
to phases, and we could actually choose arbitrary phases 
7t.4) ^ Iwta)- Choosing = |o/) Ti|ls) 

fixes those phases and actually selects the s-channel over 
the t or u channels. 

We can now define a pure spin network state which 
maps its quantum fluctuations on the thermal fluctua¬ 
tions of a given classical statistical model -for instance 
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the Ising model- by 

|V-)= ^ 

{cr„=±} 

with yl[a„] + 5 

where we sum over all the possible spin ay configurations 
modulated by an arbitrary phase ©[cti,] and an ampli¬ 
tude A[ay] which has been chosen to be an Ising nearest- 
neighbors model with a coupling constant J and mag¬ 
netic fields By. We have mapped the Ising spin tT„ = ± 
onto the space orientation of the square volume eigen¬ 
values The state '0'^’'^/^ represents a particular 

configuration of the spin network and the full state is a 
quantum superposition of them all. Defined as such, the 
state is unnormalized but its norm is easily computed to 
be the Ising partition function .Rising: 

Zl,ing = ^ e-^ ^ e-^ (8) 

{< 7 „} { ct „} 

using the usual condensed matter notation (u,u') for 
nearest neighbor vertices or the usual loop gravity nota¬ 
tion s(e),t(e) respectively for the source and target ver¬ 
tices of every (oriented) edge e of the graph. 

In principle we could define such states with any ampli¬ 
tude involving Ising spins and generalize our Ising spin 
network states to any other condensed matter models 
built from 2 -level systems. 

The intertwiner states living at each vertex are now 
entangled and carry non-trivial correlations. More pre¬ 
cisely this state exhibits Ising correlations between two 
vertex i, j: 

(cr.CTj) = —^ ^ cr,crj |A[cr„]|2 (9) 

Those correlations are between two volume operators at 
different vertex which are in fact components of the 2 - 
point function of the gravitational field. So understand¬ 
ing how those correlations can behave in a non-trivial way 
is a first step toward understanding the behavior of the 
full 2 -point gravity correlations and for instance recover 
the inverse square law of the propagator. 

We see that in order to emulate correlations of a clas¬ 
sical statistical system, the amplitudes appearing in the 
pure state must be the square root of the classical Boltz¬ 
mann factor. More generally, a spin dependent phase 
0 [ct„] could be introduced in the definition of the state 
itself as done in Eqj^ This phase could for example con¬ 
sist in complex valued local magnetic fields, which occur 
in the Lee-Yang zeroes theorem and are relevant to some 
decoherence models [20]. For instance, if we change the 
phases in the definition of the two states this would 

clearly change the Ising spin network state defined above, 
but the phases would entirely be re-absorbed into the 
Ising Hamiltonian as purely imaginary magnetic fields. 
Such a phase term doesn’t actually change the spin cor¬ 
relations, (•)e = (-(e^o- They will nevertheless affect the 


expectation values of dual operators that shift the spins 
and will affect the Hamiltonian constraints satisfied by 
the state as we will see below. For now, we consider this 
phase to be equal to unity for the sake of simplicity. 

Similar states have been studied in a quantum informa¬ 
tion approach to condensed matter physics and appear 
to have nice properties with respect to entanglement en¬ 
tropy or even for quantum computation purposes mM- 
The spirit of this approach consists in constructing quan¬ 
tum states with well controlled physical properties and 
then understand them as a ground state of a particular 
dynamic. This is the same perspective we are advocating 
for loop quantum gravity here. 

B. Low and high temperature expansion definition 

In order to better understand the structure of the Ising 
spin networks and how to build them from the basic loop 
quantum gravity operators, we look at them for the per¬ 
spective of the low and high temperature loop expansions 
of the Ising model. This is also interesting from the point 
of view that the duality between the loop expansions at 
low and high temperatures allows to characterize exactly 
the critical point of the Ising model, especially for the 2d 
square lattice, which is self-dual (see e.g. (Hj). 

a. Low temperature expansion The partition func¬ 
tion of the Ising model admits different representations 
which are more relevant in different temperature regimes. 
We restrict ourselves here to the case with no local mag¬ 
netic field. In the low temperature regime, typically be¬ 
low the critical temperature, the ground state of the sys¬ 
tem is an ordered phase. It is then natural to extract the 
ground state contribution from the partition function and 
focus explicitly on the excitation contributions. Such a 
representation is called the low temperature expansion, 
or cluster expansion. For a planar graph, as the square 
lattice, it reads 

Rising = 51 (10) 

= 2e^^ 5]] = 2e^^ 

CgC 7*65* 

where z is the valence or order of the graph’s nodes, fixed 
here at z = 4 for the square lattice. We are summing over 
all clusters C, that is all subsets of vertices of the graph 
r. The expansion is obtained by distinguishing the up 
spins from the down spins. A cluster C is equivalent to 
a even subgraph 7 * C T* of the dual lattice, that is such 
as the valence of each node of the subgraph is even, i.e. 0 
(the node does not belong to the subgraph), 2 (the sub¬ 
graph goes through that node) or 4 (which is the maximal 
value). This even subgraph on the dual graph is simply 
the boundary of the cluster, as illustrated on figure 
Pc = Pj* is the number of edges of the dual subgraph 
or equivalently the total perimeter of the cluster. One 
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can decompose an even subgraph in terms possibly in¬ 
tersecting loops, with loop intersection corresponding to 
the 4-valent nodes of the subgraph. 
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FIG. 4. An example of a cluster configuration C C F com¬ 
posed of two groups of Ising down spins amidst the up 
spins u-j-. The cluster boundary is made of two loops on the 
dual lattice, which form an even subgraph. 


This low temperature expansion leads to an alternative 
definition of the Ising spin network, which we note IV’lt)- 
We start from the totally ordered state \u^u^ ■ ■ ■) and act 
with an operator that flip all spins inside a given cluster 
C. This reads: 


c 

\C) ( 11 ) 

vec 

where f is the notation used for the Pauli matrices. Tx is 
the flip operator that sends the up spin |u^) on a down 
spin |m^) and vice versa. It can be constructed in terms 
of geometric operators, from area and scalar product op¬ 
erators acting on the intertwiner space, as shown in ap¬ 
pendix 1 ^ 

b. High temperature expansion In the high temper¬ 
ature regime, their exist various expansions and we focus 
on the loop expansion of the partition function. It is 
constructed using the following well-known linearization 
identity for Z 2 variables: 

e"^*^ = (cosh J -I-cr sinh J), Vcr = ±1. (12) 

Expanding the nearest neighbor exponentials leads to a 
loop expansion of the Ising partition function: 

■Rising — (cosh T) e)crt(e) tanh J) 

{ct„} e 

= 2^(cosh J)^ y^(tanh J)^~' . (13) 


network state IV’ht) 

E 

|'0ht) = ( cosh - I tanh 




7cr ' 


17 ) 


11^7 

where we have switched the intertwiner basis to: 


l±) = (l«t) ± \ h )) ■ 


(14) 


(15) 


We are now summing over all possible subgraphs 7 (de¬ 
fines as arbitrary subsets of edges of the initial graph 
r) and not restricting ourselves to even subgraphs. The 
index e)) actually registered the parity of the valency of 
a vertex v with respect to the subgraph: e)) = -|-1 if u 
is attached to an even number of edges of the subgraph 
7 , while e 2 = — 1 if it is attached to an odd number of 
edges. It is a non-trivial consistency check to show di¬ 
rectly that the norm (V'htIV'ht) is actually simply the 
Ising partition function. 

This high temperature state is constructed not form 
the full spin-up state but from the plus state |-|-) = 
^ (|u^) + |u^)) which actually sums over all possible 
combinations of spins up and down. This seems natu¬ 
ral since the high temperature regime of the Ising model 
is totally disordered. 

From the perspective of loop quantum gravity, the |-|-) 
state is actually |0s) and the operator acting at a ver¬ 
tex V is simply the normalized squared volume operator 
= -^Uy. This allows to define this high temperature 
expansion of the Ising spin network entirely in terms of 
geometric operators. The interested reader will And more 
detail in appendix [X| 


C. Loop representation 

Loop quantum gravity is based on the reformulation of 
general relativity as a SU(2) gauge field theory. Its basic 
observables, and then operators at the quantum level, are 
the holonomies. From the viewpoint of the spin network 
wave-functions, the gauge invariance is ensured by the 
intertwiners. It is always enlightening to understand how 
to reconstruct some spin network states from the basic 
holonomy operators and we will investigate below how 
to derive such a loop decomposition for our Ising spin 
network states. 

So the Ising spin network wave-functions for a vanish¬ 
ing phase 0 [fT„] = 0 reads: 


We are summing over all even subgraphs 7 of the initial 
lattice r. These subgraphs, as said earlier, can be seen 
as sets of possibly intersecting loops on the square lat¬ 
tice. Then P-y is the number of edges, or equivalently 
the perimeter, of 7 . We proceed similarly from 0 and 
we define the high temperature form of the Ising spin 


'^(de) = E ^ G) CTi, I (16) 

[o-„] V e V ) 

where are the spin-1/2 Wigner matrix represen¬ 

tation of the SU{2) group element g^ associated to every 
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oriented edge e. Changing the orientation of an edge 
simply simply switches to its inverse Using the 
relation between the eigenstates of the squared volume 
operator and the spin zero intertwiner states in the three 

pairing channels, + e^^ |0t)eT^ |0„)j, 

we obtain a decomposition of the Ising states over tessel¬ 
lations T on the lattice: 


n Xl/2 fn^e) 

r [cr„] Cer \eec / 

(17) 


where the phase 0(T, [cti,]) depend on both the tessella¬ 
tion and the Ising spins, and Xi /2 is the character (i.e the 
trace) of the fundamental representation. We call here 
tessellation a set of loops covering every link of the lat¬ 
tice once and only once (or equivalently a partition of the 
edge set of the lattice in terms of closed loops). Defining 
an auxiliary variable = 0,1, 2, as a value for the three 
channels s, t, u respectively, see figurej^ the phase takes a 
simple form, 0(T, [cti,]) = 6 '([t«], [cr«]) = - ^vC^v Here, 
we have omitted a factor (2/3)^ in the normalization of 
the wave-function ( |17[ ) coming from the decomposition 
of the up and down states in terms of the spin-0 

states in the s, t and u channels. 

It is interesting to note that the amplitude associated 
to a given tessellation is a partition function of a 2D Ising 
model with local imaginary magnetic fields. For a given 
tessellation, or equivalently channel values for every 
vertex, the effective Hamiltonian is 


TFioop — 2 ^ ^ ^w ^ ^ 

{v^w) 


(^By + 


(18) 


An analytical solution of the 2d Ising model with mag¬ 
netic fields is actually still unknown. In one dimen¬ 
sion, for By = ih, there exist an infinite number of cou¬ 
ples {J,h) for which the correlation length diverges |23| . 
In fact, studying models with imaginary magnetic fields 
gives information on the onset of phase transitions, relax¬ 
ation or decoherence timescales |24H26| . Indeed, the Lee- 
Yang theorem explores the zeros of the partition function 
which in turn are related to the divergence of thermody¬ 
namical quantities like the free energy. It would be in¬ 
teresting in future investigation to look at that potential 
interplay between the Lee-Yang zeroes and the structure 
of the Ising spin network states in the thermodynamical 
limit. 

Figure gives an example of a particular tessellation 
which as defined above is a set of loops passing once and 
only once through each link of the lattice. As a potential 
candidate for the generic solutions of the Hamiltonian 
constraint, this loop decomposition gives us information 
on the loop structure of physical states. Generically, the 
emergence of non trivial correlations is a consequence of 
the quantum superposition of all possible sets of extended 
loops covering the spin network. 



FIG. 5. One example of a tessellation that appears in the loop 
decomposition of the Ising spin network state (for clarity some 
loops have been omitted at the boundary). The thick loops 
are examples of extended loops contrary to the dashed ones. 
It is those that are at the origin of the non trivial correlations. 


Having the loop representation, it is interesting to 
check the results mentioned in the previous section and 
see that, apart from a normalization factor, we have 
the same expressions of the norm and correlations. The 
method here is to write all the loop contributions as a 
product at each vertex u of a matrix element {ay\M\ay). 
This matrix is quite reminiscent of the transfer matrix 
method used in statistical physics. In our case, it is pro¬ 
portional to the identity as is expected from the very 
construction of the state. 

We give here some details on this calculation. Even if 
this specific evaluation can be done in a more straightfor¬ 
ward way, the purpose here is to illustrate the method. 
We write the norm Af 


A7=E E 

r,r' {cr„.<} 

n xi/2 n n ^1/2 n 


\eGC J C'GT' 




^ ^ A*[a'y]A[ay] 


(19) 


with at = Y^y^vty, = {l^yty\ (g)« ty) = Yly{ty\ty). 

We then basically define the matrix M as follows 


n(^;|M|a„)= E (20) 

This matrix M is obtained by multiplying a phase matrix 
P and an intertwiner matrix / whose elements are simply 
the different overlaps (t(,|t„), M = *PIP, such that 


/ 1 1/2 -l/2\ 

7=1/2 1 1/2 

V-1/2 1/2 1 J 



( 21 ) 


Doing this calculation gives M = |l where we conclude 
that the norm is proportional to the Ising partition func- 




















tion (the proportionality factor being unimportant for 
the correlations we are interested about), 

( 22 ) 

The same computation gives the expected correlation 
functions for the intertwiner spins. 


D. Density 

With this loop decomposition of the Ising spin net¬ 
works at hand, we can look at the structure of its cor¬ 
responding density p{ge) = IV'(5e)Pj which defines its 
probability profile. 

Squaring the Ising wave-function '!/’(5e) dressed with 
spins 1/2 on all edges of the lattice, we get tensor prod¬ 
ucts of two spins 1/2 on each edge, which recouple to 
a spin 0 or a spin 1. Therefore the probability density 
p{ge) = IV'(5e)P will once again decompose as a gauge- 
invariant function on SU(2)^ as a superposition of spin 
networks dressed with spins 0 or 1 on the edges and ad¬ 
missible intertwiners at the vertices. 

The spin configuration consisting in only spins 0 ev¬ 
erywhere on all edges actually corresponds directly to 
the calculation of the norm N of the density computed 
above. Up to a normalization factor, we have shown that 
we recover the Ising partition function. We generalize 
this method to compute the probability of the spin con¬ 
figurations, which also involve some spins 1, and we will 
show how they are related to correlation functions of the 
Ising model. 

Using the same notation as previously, namely \t) 
which represented the tensor product of a particular con¬ 
figuration of intertwiner states depending on some given 
values ty, we can consider the probability density p{ge) 
as a spin network state itself. 

|p(5e))= E (23) 

We would like to compute the overlap between the state 
\p{ge)) and the state |5) associated to a chosen configu¬ 
ration S of spins and intertwiners: 

(5|p(5e))= E (24) 

(7,(7' ,t,t' 

This defines a transfer matrix between the Ising spins 
ay and cr/: 

{(^ya'y\M^\®ya'y)= e't("'‘'-"‘)(5|t,t') (25) 

{t.K} 

We cut the spin configuration S into little pieces Sy defin¬ 
ing the spin configurations around each vertex v. The 


various cases are given on figj^ Then we realize that the 
transfer matrix can be factorized vertex by vertex: 

((g)„cr/|M‘^| (g)y a'y) = W{(j'y\M^\ay) (26) 

V 

with (cr/IMf |cr.„) = . (27) 


Omitting the unambiguous vertex index v, we proceed 
as in the previous section and write the transfer matrix 
as =‘ PI^P where P is the phase matrix defined in 

the previous section in eqn.(21l and the 3x3 matrix 
interpolating between ty and t/ and whose elements are 
the scalar products {Sy\ty,ty). 

The task is then to evaluate those matrices. On the 
one hand, we have the states which are simply 

tensor products of the states |0, l)s and in the two other 
channels as defined in section II A[ eqn.(A2). On the 
other hand, we have the spins 0 and 1 and we need the 
intertwiners between spins 1. There are three different 
cases, whether we recouple between 2,3 or 4 spin-1 at 
the considered vertex, as illustrated on figure The 
bivalent and trivalent cases define a unique intertwiner. 
The four-valent is more involved since the intertwiner 
space is 3-dimensional. 


0 1 
0 - i - 1 0 - i - 1 


1 1 


1 1 



1 


FIG. 6. Graphical representation of the spin 1 intertwiners 
which are used in to extract parallel transport information 
from the density. 

For the purpose of clarity, we give explicit details on 
a particular example. The other cases can be treated 
similarly. We focus on the case with spin 1 on the third 
and fourth edges, as the first case of figure]^ This inter¬ 
twiner, denoted |/ooii)j is equal to | 0 i 2 )|l 34 ) , with 


IO12) 

II34) 


7l( 

1 

7! . 


t 

t i 
t i 


t 






The column notation mirrors the presence of two copies 
of spin 1/2, the top line representing the state one copy 
and the bottom one the other. The matrix /qoii then 
consists in the scalar product (0i2l34|0s_t_„0s_t_„). Cal¬ 
culating those overlaps is now straightforward. The in¬ 
terested reader will find the explicit expressions of the 
states in appendix. We obtain the intertwiner 

matrix /qoh: 


^0011 — — 


4 


V-1/2 


1/2 

0 

- 1/2 


-l/2\ 

-?) 


(28) 
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The last step is to perform the multiplication by the 
phase matrix P to finally have the specific 2x2 transfer 
matrix Moon, 

-^0011 =-^ (^2 ^ 


Let us remember that this matrix acts as an operator 
insertion between the Ising spins and a'y Looking at 
the probability amplitude {S\p{ge)) defined in eqn.(24l, 


we will have to insert all of the transfer matrices for each 
vertex. We give in appendix]^ the complete set of trans¬ 
fer matrices and operators associated to each type of con¬ 
figuration of spins and intertwiners. Having those, we can 
see that for a particular configuration S the amplitude is 
given by a correlation function of the Ising model: 

(5|p(5e))= E A*K]A[a^]Yl{a'JM^\a^), (30) 


where each of the operator insertion {a[j\M^ \ay) will con¬ 
sist in a Pauli matrix combination, switching or phasing 
the Ising spins at every vertex v. 

This shows how to express the probability density de¬ 
fined by the Ising spin network in terms of the Ising cor¬ 
relation functions. 


III. HAMILTONIAN CONSTRAINTS 
A. Definition and algebra 

Up to now, we have focused on defining Ising spin net¬ 
work states such that their n-point functions map onto 
the correlations of the standard Ising model and on how 
to generate them using geometric and holonomy opera¬ 
tors, but they remain kinematical states of loop quan¬ 
tum gravity. It would be more interesting and physi¬ 
cally relevant if we could build physical states solving 
the Hamiltonian constraint operators or if we could in¬ 
terpret our Ising spin networks as at least approximate 
physical states in some regime. Indeed, the Hamiltonian 
constraints are crucial in loop quantum gravity, they gen¬ 
erate the dynamics of the theory by implementing its in¬ 
variance under space-time diffeomorphism and are the 
quantum Einstein equations for quantum gravity. How¬ 
ever there isn’t yet a systematic method to determine 
physical states, or at least a perturbative approximation 
scheme, in loop quantum gravity despite several recent 
lines of research (especially using spinfoam models |27p 
since Thiemann’s original proposal [5S] (see e.g. pS] for 
a review of the various approaches through toy models). 

Here we will build some local Hamiltonian constraint 
operators using the basic geometric operators of loop 
quantum gravity, such that the Ising spin network states 
are their only unique solutions. Of course, these do not 
have any a priori link with actual gravity or any proposal 
of dynamics for loop quantum gravity. Moreover, we are 


working on a fixed graph and one might expect the quan¬ 
tum gravity dynamics to act on the spin network graph. 
Nevertheless, despite all these shortcomings, there might 
be some lessons to draw for such a toy construction. 

First it illustrates the type of Hamiltonian constraints 
that would lead to Ising-like spin network states and thus 
to admitting non-trivial long range correlations in a con¬ 
tinuum limit. Second we will build our constraint opera¬ 
tors from the volume operators acting on nearest neigh¬ 
bor vertices, or equivalently around loop on the dual lat¬ 
tice. This suggest a change of perspective from the usual 
methods to construct regularized Hamiltonian constraint 
operators in loop quantum gravity, which focus on the 
holonomy operators around loops on the actual lattice. 

Let us thus construct some Hamiltonian constraint op¬ 
erators characterizing the Ising spin networks. From the 
condensed matter perspective, a first approach could be 
to look at the parent Hamiltonian technique and build 
a Hamiltonian operator whose ground state is our con¬ 
sidered Ising state [5T]. We will follow a simpler path, 
looking at the action of the basic Pauli matrix operators 
on the Ising spins and using them to define simple local 
Hermitian Hamiltonian constraint operators characteriz¬ 
ing uniquely our Ising states. 

Let us consider the Pauli operator acting on the 

Ising spin at the vertex v and first focus on the switching 
operator |u^) = |u^), |m^) = |u^). Its action on the 

Ising state can be re-absorbed in the amplitude: 

TvH) = (31) 

{a-i} i^v 

= ^ A[a,] (g) |a.) 

{iTi} i 

with A[a^] = ei + 5 . 

Up to the magnetic field term, the change of sign of the 
Ising spin at the vertex v translates into the insertion of 
an additional factor involving its nearest neighbor spins, 
^ ■\Yg turn this factor in its linearized form, 
g-Jo'„cr„ _ gosh(J) — (T„tTu, sinh(J), which allows us to 
define the following Hamiltonian constraint operator at 
the vertex v. 

h: = t: - ft . n ( 32 ) 

{v^w) 

where we use for simplicity the auxiliary operators : 

ft. = cosh By — Ty sinh By (33) 

hyw = cosh J — fyT^ sinh J (34) 


^ These auxiliary operators are easily invertible: 

h~^ = cosh^t; + sinh St, 
h~^ = cosh J + sinh J 
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We have just showed that these constraints annihilate 
our Ising spin network state: 

h:\^) = 0. (35) 

Let us remember that the operator is simply the nor¬ 
malized square volume operator acting at the vertex v. 

Similarly, the action of the other Pauli matrix oper¬ 
ator leads us to define another set of Hamiltonian 
constraint operators 

Hy = tI + \tIK n (36) 

Since the operators can all be constructed from 

the area and volume operators, those Hamiltonian con¬ 
straints are entirely expressed in terms of geometric op¬ 
erators. This differs slightly from the usual construction 
of the Hamiltonian constraint operators in loop quantum 
gravity, based on Thiemann’s trick |28| and which in¬ 
volves the holonomy operator (as regularization for the 
curvature) around loops of the spin network graph. Here, 
the product of operators over the nearest 

neighbors of the vertex v can be also considered as living 
on a loop but on the dual lattice, as illustrated on fig|^ 
We will discuss more this shift of perspective below in 
section IHI Cl 


'VW2 



'VW4 


FIG. 7. A nearest neighbor vertex operator can be seen as a 
loop operator on the dual lattice around the vertex v. 

We note here that the possible local spin dependent 
phases that could be added in the definition of the state 
(Eqj^ which do not change the correlation functions 
would change the Hamiltonian constraints previously de¬ 
fined. 

By construction, the Ising spin network state is solu¬ 
tion of all these constraints, = 0. To check that 

there are not any other constraints, we check the algebra 
they generate. The constraints for two vertices which are 
not nearest neighbor vanish. The commutation relations 
on a single site v or for nearest neighbors {v, w) read: 


±±V 5 

— fry] 

' V V 

(37) 

TJX TJX 

= 2sinh J hyyyT^^Hl^ 


Ry Ry 

, -L-Lyj 

= 2sinh J hyyyT[^H!l^] 



= 2sinh J hyyyT^^H^] 



Every other commutation relations are zero. So the con¬ 
straint algebra does not generate any further constraints 
satisfied by our Ising spin network states. 

Our Hamiltonian constraints Hy exist for all values 
of the Ising coupling J. They do not depend on the 
specific phase of the Ising model, it doesn’t see a priori 
the structure of the correlations or the phase transition. 
We think that this setting would be the perfect testing 
ground for any coarse-graining scheme for spin network 
states in loop quantum gravity (e.g. [T^ [30H33])- In¬ 
deed the phase diagram, coarse-graining and phase tran¬ 
sition of the 2d Ising model is entirely under control and 
we know what to expect from the coarse-graining flow. 
The coarse-graning procedure for loop quantum gravity 
should reproduce the same flow on the correlations of the 
Ising spin network states. It will then be enlightening to 
understand what happens at the level of the Hamiltonian 
constraint operators, how the emergence of large scale 
correlations is taken into account and what triggers or 
signals the phase transition at the critical Ising coupling. 

Another approach would be to have Hamiltonian con¬ 
straints specifically tuned to the coarse-graining prop¬ 
erties of the Ising partition function, which would not 
determine the Ising spin network states for arbitrary val¬ 
ues of the coupling J but that would select specifically 
the critical point. Such Hamiltonian operator would re¬ 
flect the exact coarse-graining flow of the Ising model. 
This would be much more complicated to realize than 
our present proposal. It would either involve a graph 
changing dynamics or implement a self-duality property 
of the Ising partition function at the critical tempera¬ 
ture, probably through the square-star or triangle-star 
relations m- We believe that this is a very interesting 
line of research, but yet out of the scope of the present 
work and we postpone it to future investigation. 

B. Unicity 

So far, we have introduced a set of constraints 
Hy’y\ijj) = 0 of which the Ising spin network state is 
a solution. We now show that is it the unique solution of 
those constraints (up to a global phase factor). 

The state of the intertwiners has the following general 
form 

1^) (38) 

[o-„] V 

where Q:[o-„] are the coefficients of the state on the in¬ 
tertwiner basis. Imposing that the state solves the con¬ 
straints H!l;'y\tj}) = 0 leads us to a detailed balance type 
condition between a configuration and another with only 
one spin flip difference at a vertex v, 

acTi,... ,o-jv = .<tjv (39) 

The solution of this relation is found as follows. Begin¬ 
ning at a particular configuration, namely the one with 
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all spins up we get the amplitude of an arbitrary 
configuration by Hipping the relevant spins in that bal¬ 
ance equation. The exponential factor that appears only 
depends on the number of pairs of anti-parallele spins. 
Using the low temperature expansion notation of section 
|II B| with down spin clusters C, we obtain that 


a[o.^] = ae , 



(40) 


where the amplitude a associated to the totally ordered 
configuration with only spins is fixed by the normal¬ 
ization condition . We recognize the low temperature ex¬ 
pansion of the Ising partition function. So we conclude 
that the amplitude is then proportional to the Boltzmann 
factor of the nearest neighbors Ising model with coupling 
J/2. Thus the constraints have a unique solution up to 
a global phase given by the Ising spin network state. We 
note that because of the relation between the constraints 
and Ijy, we only require either one to reach this con¬ 
clusion. For more complex qubit models, this might not 
be the case. 

This discussion also enlightens the action of the Hamil¬ 
tonian constraints and the dynamics they could create. 
Their role is to impose a detailed balance condition be¬ 
tween two different configurations of intertwiners. Such 
behavior is often encountered when studying stochastic 
systems relaxing toward equilibrium such as for instance 
Glauber dynamics for Ising models which consists exactly 
at looking at local spin flips. 


C. Insight for Loop Quantum Gravity 

The primary purpose of our Ising spin network states 
is to provide a toy model framework to study coarse- 
graining in loop quantum gravity. Working in the con¬ 
trolled environment of the Ising model with its explic¬ 
itly known correlations seems to be a perfect testing 
ground to investigate coarse-graining procedures, contin¬ 
uum limit definitions and phases transitions. But beyond 
this aspect, it turns out that they could also bring some 
insight into the structure of the dynamics of loop quan¬ 
tum gravity. 

Indeed, as we have pointed out earlier in section [III A[ 
the Hamiltonian constraint operators that we introduce 
for the Ising states are different from the typical con¬ 
struction in loop quantum gravity, based on holonomy 
operators around loops of the spin network graph or cre¬ 
ating such loops |281 IM] (also see [3S] for a more recent 
reformulation of the holonomy and grasping operators 
in spinorial terms). These holonomies usually enter the 
constraints as a regularization of components of the cur¬ 
vature tensor. Here the natural structure of our Hamil¬ 
tonian operators involve volume operators on dual loops 
(living on the dual graph) acting on all nearest neighbors 
of a given vertex. Although our present construction has 
clearly no a priori link with gravity, it seems closer to 


the intuition of gravitational waves deforming volumes 
and shapes from vertices to vertices. This suggests to 
look for a reformulation of the loop quantum gravity dy¬ 
namics in such terms, or more generally to investigate 
the relation (e.g. under the form of a dual expansion) 
between the two types of Hamiltonian. This might help 
seeing gravitation waves (or at least geometry deforma¬ 
tion waves) emerge in loop quantum gravity directly at 
the level of the Hamiltonian constraint algebra and not 
only in a large scale semi-classical limit. 

From our viewpoint, this would require, first un¬ 
freezing the spins on the spin network edges to allow 
from holonomy operators to shift those spins, second to 
work out how our Hamiltonian constraints for the Ising 
spin network states can be re-written as some fixed-spin 
projection of more general operators built of holonomies. 
To this purpose, our decomposition of the Ising states in 
terms of basic holonomies given in section jll Cj might be 
a good starting point. Then we hope to generalize this 
discussion to the full loop quantum gravity framework 
outside the very restricted toy model of the present Ising 
spin network states. 


IV. PHASE DIAGRAMS AND CONTINUUM 
LIMIT 

A. Orientation alignment 


The effective spins we used were defined using the two 
eigenstates of the square volume operator which physi¬ 
cally correspond to two different possible orientations of 
the fundamental volume. It is interesting to identify a 
regime in which all orientations would be aligned, for in¬ 
stance either all in the positive sector of the C/„ operators 
or in their negative sector. This is exactly what happens 
for the Ising model, in both 2d and 3d, at low tempera¬ 
ture or equivalently at high enough coupling J in our set¬ 
ting. Indeed beyond the critical coupling, 
for the regular square lattice na, the Ising mode pre¬ 
dicts an ordered state, with all the Ising spins align with 
each other. More precisely, a phase transition occurs for 
the Ising model from a disordered to an ordered state. 
Figure shows the generic behavior of the magnetiza¬ 
tion (cr„) as a function of the coupling J (see eqn.(42l 
for exact formulas in the 2d case). Above the critical 
coupling, the system acquires an average orientation di¬ 
rection, randomly picked when passing the phase transi¬ 
tion. Of course the system is perfectly ordered only in 
the infinite coupling limit (zero temperature in statistical 
physics) and gets quantum fluctuations away from it. 

As for the Ising spin network states, in the infinite 
coupling limit, the state reduces to the sum over the two 
limit states, the one where the intertwiners are all in the 
positive volume orientation state plus the one where the 
intertwiners are all in the negative volume orientation 
state. As the coupling decreases, opposite spin clusters 
appear and their typical size increases with the corre- 
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FIG. 8. Generic representation of the behavior of the mag¬ 
netization (in normalized units) as a function of the coupling 
passing through a second order phase transition. 


lation lengths. Towards the phase transition, we have 
clusters of all sizes with no apparent ordering or pre¬ 
ferred orientation. Beyond the phase transition, at lower 
coupling, we are in the disordered phase, with no ori¬ 
entation alignment. An ordered universe, with aligned 
volume orientation, would therefore live in the ordered 
phase at high coupling J (i.e low temperature). 


B. Distance from correlations 


In a non-perturbative background-independent ap¬ 
proach to quantum gravity, the very notions of distance, 
locality or metric need to be reconstructed from scratch 
in the absence of a background geometry. Understand¬ 
ing the emergent geometry from the quantum state is 
still an ongoing issue and we hope to be able to recover 
all this geometric information from the correlations and 
entanglement in the quantum state m- 

So the idea we pursue here is that a notion of distance 
emerges from the correlations in the spin network state. 
Using the known results on the Ising model, we get the 
Ising spin network correlations in terms of the natural 
lattice distance and then work backwards attempting to 
see to what extent the distance between two spins can 
get be extracted from the sole knowledge of the corre¬ 
lations between those spins. By construction, the Ising 
spin networks have exactly the same known correlation 
behavior than the classical Ising model, namely 


(cTjCTj) - {(Ji){(Tj) OC 


Cexp^J^, 


J^Jc 

J=Jc 


(41) 


where ^ is the correlation length (in units of the lattice 
spacing) and C a positive constant. Here the distance 
between the vertices \i — j\ is the graph distance as natu¬ 
rally defined on the 2d regular square lattice. The exact 
form of the magnetization and the behavior of the corre¬ 
lation length near the transition are known but won’t be 


used in the following]^ 

((TiCTj) OC |j_^-|i/4 

vX * 

Ordered Disordered 
PHASE PHASE 

0 

FIG. 9. Phase diagram with renormalization flow of the Ising 
spin network state with the algebraic decaying correlation 
function at the critical coupling. 

It is possible to invert the correlation formula so as to 
express the distance as a function of the correlation itself. 
The distance will be generically a monotonous decreas¬ 
ing function of the proper correlation between the two 
spins, {aiaj)p = {piOj) — {ui){uj)^ but the precise func¬ 
tion will depend on the considered phase and regime. 
Away from the critical point, we have an exponential de¬ 
cay of the correlation and we can define the distance as 
d{i,j) = —^ \og{ai(jj)p up to an additive constant. The 
correlation scale ^ becomes the new length unit. At the 
critical point, things become especially interesting. The 
correlation length blows to infinity and we have an alge¬ 
braic decay of correlation. We can now define the dis¬ 
tance as d{i,j) = lf{aiaj)p. 

From the point of view of gravity, the most interest¬ 
ing case is clearly an algebraic decay of the correlations. 
Ideally, we’d like to derive a quadratic decay of correla¬ 
tion, just as in standard quantum field theory, in order 
to retrieve Newton’s gravity law. In order to get this, we 
expect of course to have to change our Ising model state 
to another better suited statistical model, but also move 
to a 3d lattice or graph structures, unfreeze the spins of 
the spin network and so on. Nevertheless, the present 
Ising spin network state allows to illustrate a couple of 
important points: 

• Reconstructing the distance from the (2-point) cor¬ 
relation depends not only of the considered statis¬ 
tical model on the specific phase of that model. 

• We naturally have a algebraic decay of the corre¬ 
lation in terms of the distance, and vice-versa, at 
the phase transition. At that point, the state also 
admits a non-trivial continuum limit. 


For reference, below the critical coupling 


{Uv)^ — 
i OC 


1 - 


1 — (tanh J)^ 


2m1/4 


2(tanh J)^ 


1 1 
J “ X 


( 42 ) 
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Finally, in order to truly validate this notion of recon¬ 
structed distance, we need in general to check the triangle 
inequality satisfied by the distances between three points, 
which would involve the 3-point fluctuations. Here we 
know that the reconstructed distance is the initial graph 
distance, which satisfies this requirement. But if we 
would consider another spin network state, not necessar¬ 
ily related to a known local statistical model, this would 
have to be checked. 


C. Coarse-graining 



One of the main challenge of the loop quantum gravity 
program is to define and understand its renormalization 
flow, from the Planck scale to large scales and derive the 
low energy behavior as semi-classical general relativity 
with a systematic method to compute the perturbative 
quantum corrections. To this purpose, coarse graining 
procedures for the spin network states and the dynamics 
of the theory are essential tools |15l l3fil - l^ . 

Understanding coarse-graining in loop quantum grav¬ 
ity is not an easy task because of the absence of a back¬ 
ground geometry or structure and the complicated nature 
of the Hamiltonian constraints. In the present frame¬ 
work of the Ising spin network state, we have fixed a 
background graph, to a regular 2d square lattice up to 
now (we will deal with its three-dimensional generaliza¬ 
tion below in the next section IVD I, and we can use this 
structure to define the coarse-graining flow as in standard 
statistical physics and condensed matter models. Thus 
this provides a neat toy model to test all coarse-graining 
procedures of loop quantum gravity. 


In statistical and condensed matter physics, renormal¬ 
ization group and coarse graining methods |36| are widely 
used to understand critical phenomena. We wish to im¬ 
port these methods and results to quantum gravity to 
understand better the emergence of critical regimes. The 
usual formulation of loop quantum gravity focuses on the 
Hamiltonian constraints, their algebra and their flow gen¬ 
erating the time evolution. Here, for our Ising spin net¬ 
work states, the Hamiltonian constraints that we intro¬ 
duced do not see a priori the phase transition: nothing 
obvious happens at the level of the algebra of the con¬ 
straints when the coupling reaches its critical value. The 
information about the critical regime and phase transi¬ 
tion is truly in the coarse-graining flow. We could, for 
example, apply the tensor network renormalization tools 
to our Hamiltonian constraints, such as it was done for 
the so-called spin-net toy models in m, in order to derive 
their coarse-graining flow and check if it fits as expected 
the known flow of the 2d Ising model. 


Since the spin correlations of the Ising spin network 
states are the same as the Ising model, we know their 
coarse-graining flow from the standard techniques in sta¬ 
tistical physics. For instance. Figure [TO] show its renor¬ 
malization group flow using the basic isotropic decima¬ 
tion method: we plot the effective coupling J' at larger 


FIG. 10. Renormalization group flow for the 2d Ising model 
using the method of decimation. The critical point is totally 
repulsive. 


scale (including the nearest and next nearest neighbor 
interactions) as a function of the smaller scale coupling 
J. We can follow the renormalization flow of the effec¬ 
tive coupling by iterating this map and the critical point 
is identified as its non-trivial fixed point. In the criti¬ 
cal regime, the state presents scale invariance which ren¬ 
ders irrelevant the micro structure of the theory and al¬ 
lows for a well-defined continuum limit. However, this 
fixed point is repulsive for the Ising model: as we coarse- 
grain, the effective coupling runs aways from the crit¬ 
ical regime towards either the low or high temperature 
fixed points. We hope to later identify a better statistical 
physics model such that the coarse-graining flow would 
run instead under the coarse-graining flow towards the 
critical regime with algebraic decaying correlation and 
continuum limit. The corresponding spin network state, 
built along the same lines as our Ising spin networks, 
would then be better suited to the definition of the con¬ 
tinuum limit and semi-classical regime of loop quantum 
gravity. 


D. Going tridimensional 


So far, we restrain the discussion of the gravitational 
state to a two dimensional lattice for exposition sim¬ 
plicity but it happens that the generalization to 3d is 
straightforward. Keeping the requirement that the lat¬ 
tice be 4-valent (so as to keep a two-dimensional inter¬ 
twiner space and the map to effective qubits), the nat¬ 
ural regular lattice is the diamond lattice as illustrated 
on FiglTT] Under the usual geometrical interpretation of 
loop quantum gravity, this lattice can be seen as dual 
to a triangulation of the 3d space in terms of tetrahedra 
dual to each vertex. This can be seen as an extension of 
the more used cubic lattice better suited to loop quan¬ 
tum gravity. The Ising spin network state and the whole 
set of results which followed are then identical : the wave 


function eqn.(17l, the Hamiltonian constraints and their 
algebra eqn.(37l are the same. 


In 3d, the Ising model also exhibits a phase transition. 
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FIG. 11. Elementary cell of the diamond lattice which is the 
natural one in 3d in loop quantum gravity for a 4-valent spin 
network. 

Even if the exact solution is still unknown and is the sub¬ 
ject of active research, much can be learn of its behavior 
around the critical coupling. Like its 2d counterpart, the 
Ising spin network state has an ordered phase which cor¬ 
responds to an orientation alignment of the elementary 
chunk of space. Concerning the 2-points correlation func¬ 
tions, information about its long distance behavior at the 
phase transition or near it can be obtained using meth¬ 
ods of quantum field theory. In d dimensions we have 

m 

where K{r) are modified Bessel functions and ^ is the 
correlation length. We mention that the correlation does 
not depend in the long distance regime on the "magnetic 
field" B (even local ones Bi if we choose to define the 
state in this way). For the three dimensional case, we 
have the simple and exact expression 

((T*crj) oc —(44) 
47r|*-j| 

Figure represents the phase diagram of the 3d Ising 
spin network (in the presence of a "magnetic field"). At 
the critical coupling, the correlations have a one over the 
spacing power law decay at long distance, still not the 
inverse square law but less exotic then the 2d behavior. 
Understanding distance from correlations in 3d is much 
more vital than in the 2d setting since in 3d their ex¬ 
ist loops that can wind around a vertex without visiting 
other distant vertex (in a relational sense). 

The expression for the correlation is obtained in a mean 
field setting where only a correlation length appears. 
However statistical systems possess another length scale, 
the lattice spacing, which leads to correction of this ex¬ 
pression through the anomalous dimension. In quantum 
gravity the Planck length plays this role and it would be 


interesting to understand its influence on the behavior of 
the correlations. 

(o-iCTj) oc 

1/Jc !/•/ 

Ordered Disordered 

PHASE PHASE 

0 

FIG. 12. The 3d phase diagram of the spin network state on 
the honeycomb lattice has a phase transition whose correla¬ 
tions behave with a one over the separation law. 


V. CONCLUSION 

In this paper, we have introduced a class of spin net¬ 
work states for loop quantum gravity on 4-valent graph. 
Such 4-valent graph allows for a natural geometrical in¬ 
terpretation in terms of quantum tetrahedra glued to¬ 
gether into a 3d triangulation of space, but it also allows 
us to map the degrees of freedom of those states to ef¬ 
fective qubits. Then we were able to define spin network 
states corresponding to known statistical spin models, 
such as the Ising model, so that the correlations living 
on the spin network are exactly the same as those mod¬ 
els. For instance, we defined some Ising spin network 
states, both in 2d and 3d, on 4-valent regular lattice. In 
2d, we use the usual square lattice, while the 3d case used 
the diamond lattice. 

We want to propose those Ising spin network states 
and their generalizations (for instance to the XY model, 
Potts models and others) as toy models to test coarse- 
graining methods in loop quantum gravity and allow for 
explicit discussions on phase diagrams and transitions. 
As a start, we used them to illustrate the relation be¬ 
tween correlations (2-point functions to be precise) and 
distances. Indeed, in a background independent frame¬ 
work with no a priori geometrical structure, all the geom¬ 
etry needs to be reconstructed from the quantum grav¬ 
itational state and can only emerge from its correlation 
and entanglement structure. Here, using the exact phase 
diagram and correlation functions of the Ising model, we 
showed how the distance could be defined a posteriori 
as a monotonous decreasing function of the correlation, 
following the intuition that two systems are close if they 
interact much and consequently are strongly correlated. 
Now the explicit map between correlation and distance 
actually depend on the considered phase and regime: in 
the ordered or disordering phase, the distance would be 
minus the logarithm of the correlation, while we get an 
algebraic decay at the critical point. More generally, 
we hope that it will be possible to export the statist!- 
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cal physics methods to study the detail of the coarse- 
graining and renormalization flow of those Ising spin net¬ 
work states. 

We further described the Ising spin network states in 
details from the point of loop quantum gravity, defin¬ 
ing them in terms of geometrical and holonomy opera¬ 
tors, showing for instance that getting long-range cor¬ 
relations required holonomy operators around arbitrarily 
large loops on the graph (and not only local holonomies). 
We also introduced Hamiltonian constraints that charac¬ 
terize uniquely our Ising states. This opens the possi¬ 
bility of, on the one hand, investigate a potential refor¬ 
mulation of the actual loop gravity constraints in similar 
terms, and on the other hand, analyze the coarse-graining 
flow of such Hamiltonian constraints and algebra in the 
thoroughly-studied environment of the Ising model where 
everything is under control. 

In fact, we view our construction as a first step towards 
studying the coarse-graining of spin network states and 
the loop quantum gravity dynamics from the condensed 
matter perspective. A similar angle of attack has been 
developed also in m towards the construction of spin 
network with controlled entanglement and in m with 
the application to quantum gravity toy models of sta¬ 
tistical coarse-graining methods such as tensor network 
renormalization. We see a few lines of development of 
the present work. First, we should deepen the analy¬ 
sis of the proposal of the distance reconstruction from 
the correlations and go beyond the 2-point correlation 
to the 3-point and many-point correlation to check tri¬ 
angular inequalities and more generally to what extent 
we get a metric structure. Second, we should study the 
coarse-graining and renormalization flow of the Hamilto¬ 
nian constraints and compare it to the standard flow of 
the Ising model. Third, we should generalize our con¬ 
struction to more generic states of loop quantum gravity, 
allowing for more freedom in the spins and intertwiners, 
and therefore extending our definition of spin network 
states to more complex statistical physics models with 
more interesting critical regimes. To conclude, we believe 
that our construction of spin network states as qubit sta¬ 
tistical models offers a perfect arena as a toy model to 
investigate the structure of correlation and entanglement 
in spin networks and the emergence of phase transitions 
in loop quantum gravity, and we hope to use this tool to 
investigate further and test renormalization methods of 
the Hamiltonian constraint operators in that context. 


Appendix A: On 4-valent intertwiners between 
spin-1/2 


In the section we present a brief summary of the dif¬ 
ferent states used in the core of the discussion and their 
relations. We consider the Hilbert space V of the spin 
1/2 representation of SU(2). The natural basis of this 
two-dimensional space is denoted We now look at 

the intertwiner space between four spins, i.e the SU(2)- 


invariant subspace in the tensor product space This 
is again a two-dimensional space, which we map onto the 
Hilbert space of an effective qubit: 

dimlnvsu( 2 ) [V®"*] = 2. (Al) 

We construct a natural basis of this intertwiner space 
by recoupling the spins by pair into an internal spin, 
which can take the value 0 or 1. They are three pos¬ 
sible pairing between the four spins, which correspond to 
three different channels, which we called s, t, u in the 
main text and which are illustrated in flg|^ Recoupling 
the spins 1 and 2 together corresponds to the s channel 
and gives the following intertwiner basis states: 

|0s) = ^ (inn) + inn) - inn) - inn)) 

|i.) = ^(ittn) + inn) 

— inn) + inn) + inn) + inn) ^ ^^2^ 


We can similarly express the intertwiner states of spin 
0 and 1 in the two other channels t and u in terms of 
the It, t) states. Then we compute the change of basis 
between the three channels: 


|0t) 

| 1 *) 

| 0 „) 

in) 


1/1 /| 0 ,) 

21,73 -i) Mn) 

1 / -1 V3\ ^|0,) 


2 V-Ts -I Vin) 


(A3) 

(A4) 


The Casimir operators, given by the scalar prod¬ 
uct between the SU(2) generators, are easily calcu¬ 
lated in these various basis using the formula J.K = 

and the transfer matrices be- 




tween basis. We compute in the (|0)]^2, |l)i 2 ) basis: 


J 1 .J 2 — 


1 

4 



(A5) 


J 1./3 


1 / 0 -y3\ 

4 2 ) 


(A6) 


This allows to compute the action of the square volume 
operator {7 as a commutator of those Casimir operators: 


C = -i[Ji.J2,Ji.J3] 

.^3 / 0 1 

“'4 (^-1 0 


(A7) 


We define the effective qubit states from the eigenvec¬ 
tors of this square volume operator: 


IWd) = l'“Td) 


(AS) 
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^ (A9) 

= (|0*)±*|1,)) 

= ^eT^(|0„)T*|lu)) 

= I (|0s) +eT^|0t)eT^|0„)) 

These will be the up and down spin states for the Ising 
model, now corresponding to the two geometrical orien¬ 
tations (of the 3d space ~ K^) leading to opposite 

signs of the squared volume. 

In this qubit basis we act with the Pauli matri¬ 

ces T^’y’^ generating the su(2) algebra. The operator 
is naturally defined as the normalized square volume op¬ 
erator: 

t "" l^td) = ± l^td) = ^ ^ ■ (AlO) 

The operator is defined from the scalar product oper¬ 
ators: 

= |u^,t) = - ( 2 J 1 .J 2 + 0 \u^,i) ■ (All) 

When discussing the low and high temperature defini¬ 
tions of the Ising spin network states, we needed the 
eigenstates of r^, which we wrote |±) = (|tt^) ± |u^))/-\/2. 
In light of the definition above, these are simply the basis 
states in the s channel: 

l+) = |0.) |-) = -i|U (A12) 


Appendix B: Matrix representation of intertwiners 

From the holonomy decomposition of the Ising spin 
network state, the method developed to evaluate corre¬ 
lation functions was based on calculating product of ma¬ 
trices associated to each type of intertwiners, see section 
III Cl for an illustration of the calculation method. For 
mathematical completeness, we give the intertwiner ma¬ 
trices for those composed of two (see Figure or three 
spin one : 


The evaluation of the recoupling of three spin-one is 
rather simple: there is a unique intertwiner given by the 
Clebsh-Gordan coefficient. If we see a spin-one as vec¬ 
tor, then the invariant subspace of the tensor product of 
three vectors is given by their scalar triple product (or 
determinant), which is defined by the Levi Civita tensor 
e (totally antisymmetric tensor). The intertwiner state 
I Jill) then reads: 


= 4 E (B2) 


where the states |m) of the spin-1 representation are de¬ 
composed as tensor products of two spin 1 /2 as: 


| 1 ) 

| 0 ) 

1 - 1 ) 



(B3) 

(B4) 

(B5) 


Obtaining the matrix is then just a matter of calcu- 
lationg the different scalar products. We then have the 
associated M matrices as defined in the main text: 


=-4^ (4 - (I) (I) 

Mu = -^ Q l-kcos(^y^fy-ksin(^y^fj,^ 

Mom = - ^ ( t ) ( y ) *^4 


Is 


4 


V-1/2 


1/2 

0 

- 1/2 


l/2\ 
1/2 
0; 


it= 


4 


/ 0 1/2 l/2\ 
1/2 1 1/2 
\l/2 1/2 0 J 


In 


4 


/ 0 -1/2 -l/2\ 

-1/2 0 1/2 
V-1/2 1/2 1 J 


/m =- 


ys 1 

4 y/2 


0 1 1 
1 0 1 
1 1 0 


(Bl) 
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